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ABSTRACT 


The study of weapon delivery probabilities has historieally been focused around 
analytical solutions and approximations for weapon delivery accuracy and effectiveness 
calculations. With the relatively recent increase in modern computing power many of the 
historical expressions can be simulated quickly with similar or more accurate results than 
the historical expressions and approximations. 

In this thesis simulation methods are used to evaluate weapon delivery probability 
parameters including circular error probable, range and deflection error probable, and 
weapon effectiveness in the single and salvo weapon scenarios. Comparisons of the 
simulation results and corresponding historical practices are made to validate simulation 
techniques. 

Additionally, standard deviations in the range and deflection direction are 
extracted from weapon impact data. Using these extracted standard deviations weapon 
effectiveness, calculations are performed. 
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I. INTRODUCTION TO WEAPONEERING CONCEPTS 


In general terms, weaponeering is the proeess of determining the quantity of a 
speeific type of weapon required to aehieve a specific level of target damage, considering 
target vulnerability, weapon effects, munition delivery errors, damage criteria, probability 
of kill, weapon reliability, etc. [1] This thesis will focus on munition delivery error 
statistics and probability of hit for various scenarios. These topics are inherently random 
in nature requiring a statistical approach for analysis. The weaponeering concepts 
discussed require a general understanding of some basic statistical definitions and 
methods. This chapter provides the necessary statistical background to follow the 
analysis in the following chapters. 

A. PROBABILITY DENSITY FUNCTION 

The probability density function (PDF) describes the likelihood that a random 
event will result in a certain value contained within a defined population. Flipping a coin 
is a classic example. There is an equal probability that the coins will show heads or tails. 
Because the probability density function must account for all possible outcomes the total 
sum of all possibilities must be one. It is intuitively obvious that the probability of 
“heads” resulting from a coin toss is 1/2. The same is true for “tails”. This yields a sum 
of all probabilities equal to unity as expected. This is an example of a discrete PDF. 
Discrete meaning the data set consists of fixed values with discontinuous jumps for the 
results. A Coin flip or roll of a die are clear examples of discrete random processes. The 
results from these events can only be; 

• Coin Toss: heads or tails 

• Die roll: 1, 2, 3, 4, 5, or 6 

Of more relevance to weaponeering considerations is the continuous PDF. A 
probability can be obtained for any result within the bounds of the population. Take for 
example an aircraft dropping unguided bombs on a target. Suppose the aircraft drops 250 
bombs with impacts as shown in Figure 1. 
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2S0 Bomb Impacts 



-2100 -150 -100 -50 0 50 100 150 200 

Deflection 

Figure 1 Sample Impaet Points for 250 Bombs 


Range is defined as the direetion the aircraft is heading when the bomb is 
released. Deflection is perpendicular to the range direction with the origin of the system 
defined as the desired mean point of impact (DMPI). The impacts clearly display the 
randomness of this delivery. Two main parameters used to describe the dataset are the 
mean and variance. The mean, orx, is the average of the parameters in the dataset and 
provides a relative location of the dataset to the target. Here the mean values for range 
and deflection are displayed by the mean point of impact (MIP). The variance is defined 
as: 


1 ;=1 


The standard deviation is defined as the square root of the variance and is labeled 
a. Standard deviation represents a key parameter for characterizing weapon delivery 
accuracy and effectiveness and will be discussed extensively. 
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The randomness of the deliveries ean be evaluated further by splitting eaeh axis 
into bins and eounting how many bombs are eontained in eaeh bin. For example, for the 
impaet data in Figure 1.1 the defleetion axis ean be split into 15 ft. bins. Each bin is then 
evaluated to determine how many impacts it contains. Figure 2 displays these results. 


Impact Histogram(15 1t range bins) 



Figure 2 Impact Histogram 


This data can also be used to calculate the probability that a bomb will land within 
a certain cross-range. For example, the bin from 37.5 to 52.5 ft contains 25 impacts. 
Therefore, knowing that 250 bombs were dropped, the probability that a bomb dropped 
will be in the range of 37.5 to 52.5 is 25/250 or 10%. This technique is rather 
cumbersome however, so the histogram is replaced by a continuous expression known as 
a probability density function (PDF). [1] 


3 















1 . 


Univariate Normal Distribution 


a. Univariate Normal PDF 


There are many different PDF’s for different types of systems. The 
histogram in Figure 2 above resembles a eommon bell eurve. This bell shaped eurve is 
known as a normal distribution and represents many physieal systems ineluding many 
aspeets of weapon delivery aeeuraey. As the number of samples is inereased it ean be 
shown that the analysis of the histogram above beeomes a better approximation of the 
eontinuous expression for the univariate normal PDF given by: 


f{x) = 


1 


-{x-nf 


2a^ 


( 2 ) 


crV^ 

In Eq. (2), p is the mean value of x and a is the standard deviation. 
Approximately 68% of the datapoints for a normal distribution will be eontained within 
±a. [1] 


b. Univariate Normal Cumulative Density Function (CDF) 


The univariate normal eumulative density funetion is defined as: 


The CDF ean be thought of as the probability that a given sample will he 
in the range from -oo to some value X. Beeause this integral eannot be integrated 
analytieally table referenees are eommonly used. Standardization of the Eq (3) is used to 
allow for one table to be used for the various eombinations of p and a. The 
transformation variable z in Eq (4) is used. 



X-JU 
(j 


( 4 ) 
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Following substitution into Eq (3) and understanding that the standardized PDF has a 
zero mean and a standard deviation of 1 yields Eq (5). [1] 



( 5 ) 


If tables are unavailable, integration of Eq (3) can also be performed using 
MATLAB as shown in Table 1. Eor example, return to the previous example for the 
probability that a sample lies within the range of 37.5 to 52.5. Elsing Eq. (2) for the PDF 
and the MATEAB Symbolic Toolbox the value can be evaluated numerically. A kn own 
value of cr = 50 will be used for this evaluation. The random data generated for the 250 
bomb impact sample was also based oner = 50 . 

Table 1 MATEAB Code Eor Integrating the Non-standard Normal PDF 

% 

MATLAB CDF CODE: 

% Input Mean 
mu=0 

% Input sigma 
s=50 

% Input Start of Bin 
a=37.5 

% Input End of Bin 
b=52.5 

% Input Total Number of Bombs 
n=250 

% Creates Symbolie univariate normal PDF 
syms X 

normpdf= 1 /(s* sqrt(2 *pi)) * exp(-((x-mu). ''2)/(2* s''2)) 

% Numerieally Caleulates PDF Integral between speeified values a and b 
Prob=int(non'npdf,a,b) 

Prob=double(Prob) 

numBombs=Prob*n 

% - 

The above code results in a probability of 0.0798 that a given impact is 
contained within the range 37.5 to 52.5. Multiplying this value by the total number of 
bombs dropped results in the number of bombs contained in this range bin equal to 
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approximately 20 bombs. These values are elose to the previously ealeulated probability 
of 0.1 with 25 bombs. If the number of bombs used to ereate the histogram is inereased 
the histogram derived value will approaeh the numerieally ealeulated 0.08. For example, 
with 50,000 bombs the histogram looks like Figure 3: 


Impact Histogram (15 1t range bins) 



Figure 3 50,000 Samples Impaet Histogram 


Repeating the above MATLAB analysis with 50,000 bombs results in the 
same probability of 0.0798 (this is not dependent on number of samples) and 3,988 
samples in the 37.5 to 52.5 range. The value displayed on the histogram of randomly 
generated samples eontained in the defleetion bin from 37.5 to 52.5 is 4,004 impaets. 
Therefore, the histogram value of 4004/50000 results in a probability of 0.0801. This is 
very good agreement and demonstrates the importanee of eonvergenee when modeling 
statistieal data. Convergenee is diseussed more in Chapter II. 
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Bivariate Normal Distribution 


a. Bivariate Normal PDF 


The bivariate normal distribution ean be thought of as the eombination of 
two independent univariate normal distributions. While the univariate distribution will 
provide the probability that a bomb may fall within a eertain one dimensional bin, the 
bivariate normal distribution ean be used to determine the probability that a bomb will 
fall within a eertain area. Take the bombing example and imaging that both defleetion 
and range univariate normal distributions are used to determine the probability that a 
weapon will fall within a eertain defleetion and range distanee from the DMPI. The 
bivariate normal PDF is shown below in Eq. (6) 


y) 


1 


2o-y 2o-/ 


( 6 ) 


b. Bivariate Normal CDF 


The bivariate normal CDF is show below in Eq. (7) 


F{X,Y) 


x=X y=Y 

I I 






1 

— p 

20-/ 

2a7 


271(7 (J 

-00 y=-co X y 


dxdy 


3. Circular Normal and Rayleigh Distributions 


(V) 


a. Circular Normal PDF 


The eireular normal distribution is a bivariate normal distribution with 
zero means and equal standard deviations for x and y. This eauses Eq. (6) to reduee to 
Eq. (8) 


fix y) = 


1 

2710 ^ 


X 




( 8 ) 
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b. Rayleigh PDF 


The Rayleigh PDF is defined as the distribution of the value r defined as 


2 2,2 

r =x +y 


After some manipulation, this results in the Rayleigh PDF shown in Eq. (10) 




2<j‘- 


(J 


(9) 


( 10 ) 


c. Rayleigh CDF 


The Rayleigh PDF ean be integrated analytieally resulting in the Rayleigh 
CDF in Eq. (11) 


F{R) = \-e' 


2a^ 


( 11 ) 


B. ERROR TYPES 


There are two main types of error to be defined when diseussing unguided 
weapon deliveries. The error ean be broken down into ballistie dispersion error and 
aiming error. 

1. Ballistic Dispersion 

Ballistie dispersion error is defined as the error in the weapon delivery eaused by 
physieal ineonsisteneies between individual weapons (weight, eenter of gravity, fm 
shape/angle bias, surfaee deviations, ete.). These ineonsisteneies eause eaeh weapon’s 
ballistie trajeetory to be slightly different. The random physieal ineonsisteneies typieally 
result in weapon delivery behavior that ean be represented by a normal distribution. [1] 

2. Aiming Error 

Aiming error is the differenee between the aetual target loeation and the weapon 
system aim point. This error is also eonsidered to be normally distributed. 
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C. ACCURACY 

Numerous parameters ean be used to eharaeterize the aeeuraey of a weapon 
system. Some of the most useful parameters are eireular error probable (CEP), range 
error probable (REP), and defleetion error probable (DEP). 

1. Circular Error Probable (CEP) 

The CEP is defined as the radius of a eirele (centered on the DMPI) that contains 
50% of the bomb impacts (see Eigure 4). 



Eigure 4 Circular Error Probable 

2. Range Error Probable (REP) and Deflection Error Probable (DEP) 

The REP is defined as a length of half of a range-bin centered at the DMPI that 
contains half of the impacts in the range direction. The DEP is defined as a length of half 
of a range-bin centered at the DMPI that contains half of the impacts in the deflection 
direction (see Eigure 5). 
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3. 


Relationship of CEP to REP/DEP 


The following equations deseribe the relationship between REP/DEP and standard 
deviation in the range and defleetion direetions. 

REP = 0.6745cr^ (12) 

DEP = 0.6745cr^ (13) 

It is also observed that if the data distribution has a zero mean and is assumed to 
be eireular (standard deviations in the range and defleetion direetion are equal) the 
following relationships between a, CEP, and REP/DEP hold. [1] 

CEP = 1.1774cr (14) 

CEP = 1.7456PEP = 1.7456DEP (15) 

D. COMBINING ERROR TYPES 

Onee the statisties of the aeeuraey of a weapon are understood it is then important 
to define the probability that a given weapon, or group of similar weapons, will damage a 
speeifie target. This is the study of weapon effeetiveness and is a funetion of the weapon, 
target, and seenario. This analysis has the potential to result in a very eomplex 
ealeulation. Eor the purposes of diseussions herein it will be assumed that the weapon 
aeouracy statisties (a, CEP, REP/DEP) for the seenario are provided and the weapon area 
of effeetiveness is also provided based on known weapon/target/seenario oharaeteristies. 
An individual weapon is defined to have hit the target if its impaet is elose enough for the 
area of effeetiveness to enelose the target. Conversely, the area of effeetiveness ean be 
eentered on the target and a hit ean be defined as a weapon impaeting inside the area of 
effeetiveness. Both of these hit definitions are equivalent representations. Having these 
parameters will allow for eomparison of analytieal approximations and MATLAB 
simulation results for two speeifie seenarios: single round per oeeasion and multiple 
round salvos per oeeasion. An oeeasion is defined as an event for whieh aiming error is 
eonsidered eonstant. Eor example, an aireraft delivering two bombs (at the same target) 
on one pass is one oeeasion. An example of two oeeasions is an aireraft delivering two 
bombs on two passes, one bomb per pass. 
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1 . 


Single Round Scenario 


As discussed in section l.B there are two main error types of interest: aiming 
error (cr^,^) and ballistic dispersion (cr^^). Take any given round and assume these two 

errors are known. To properly simulate the weapon impact location both errors need to 
be accounted for as shown in Figure 6. 



Figure 6 Single Round: Aiming Error and Ballistic Dispersion 

This is an example of one occasion of a single round delivery. A four occasion 
single round delivery can be seen in Figure 7. 
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Figure 7 Single Round Seenario: Four Oeeasions 

For a typieal seenario, as shown in Figure 6 and Figure 7, it is eommon that 
and are not equal. Usually, the aiming error is a more signifieant error eontribution 
than the ballistie dispersion. 

2. Salvo Scenario 

A salvo, as one might expeet, is defined as multiple rounds per oeeasion. The 
salvo is used to inerease the probability of damage to the target. An example of a four 
oeeasion seenario with 5 bombs per salvo is shown in Figure 8. 
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Figure 8 Salvo Scenario: Four Occasions; Five Bombs/Salvo 

It is interesting to note that these two scenarios have significantly different 
calculations to determine the probability of hit. For the single round scenario the two 
error types, ballistic dispersion and aiming, can be considered independent for each 
weapon. For the salvo scenario this is clearly not the case as the aiming error for one 
salvo biases the results equally for all bombs of a given salvo. This results in a more 
complicated analysis to properly calculate the probability of at least one hit on the target 
given a salvo scenario. Examples of two approximations to salvo effectiveness are 
shown in Eq. (16) and Eq. (17). 
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Where: 


Pn(n) = Probability at least one round hits the target 
n = Number of Rounds/Salvo 

Rt= Target Radius (also can be thought of as the effective weapon radius) 
cr^ = Aiming Error Standard Deviation 


= Ballistic Dispersion Standard Deviation 

c = Adjustment Factor (Typically 0.9 to 1.0) [2] 

Monte-Carlo simulations as outlined in the following section can also be used to 
evaluate the salvo scenario effectiveness. Chapter 111 will deal exclusively with accuracy 
and probability of hit calculations and simulations for both the single round and salvo 
scenarios. 


3. Simulation Implementation 

The Monte-Carlo simulation used to generate salvo scenario results are performed 
as shown in Figure 9. 
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Figure 9 Salvo Effectiveness Simulation Flowchart [From 3] 


The process shown in Figure 9 can be easily implemented in MATFAB allowing 
for potentially more accurate results then the analytical approximations to be discussed in 
Chapter III. 
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II. TEST DATA CHARACTERISTICS 


When comparing statistical datasets with analytical solutions/approximations it is 
critical to monitor convergence. Convergence, for the purposes herein, will be handled 
using the number of samples defined for a given simulation. The simulation result is 
compared to an appropriate analytical expression if available. By studying the 
comparison of multiple simulation runs and analytical results convergence can be 
observed. If the dataset is not suitably converged the number of samples is increased and 
the simulation is repeated. This process is continued until suitable convergence is 
achieved. 

A. UNIVARIATE NORMAL 

The REP of a univariate normal distribution is shown in Eq. (18). 

REE = 0.6745cr^ (18) 


This equation is the analytical representation of the CEP. To properly compare 
this analytical representation to simulation results it is important to monitor the 

convergence of the data set. Eor example, a given dataset has a known ^ ~ which 
will result in a REP of 33.725 ft. Using a Monte-Carlo simulation to generate a normal 
dataset and extract the REP for various numbers of samples results in Table 2. 

Table 2 REP Convergence Data 


Number of Bombs 
Simulated 

REP (ft) 

Run 1 

Run 2 

Run 3 

10 

64.500 

38.363 

42.723 

lo^ 

33.105 

28.292 

26.624 

lo^ 

31.955 

33.616 

35.027 

10" 

33.369 

33.794 

33.792 

lo^ 

33.655 

33.535 

33.430 

10^ 

33.716 

33.732 

33.706 

10' 

33.742 

33.720 

33.725 
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It can be seen that even with 107 samples the simulation has not eonverged to the 
exact value for REP ealeulated from Eq. (18). However, the praetiee of weaponeering 
rarely demands this level of preeision allowing for an aeeeptable number of samples to be 
used to generate appropriately eonverged Monte-Carlo simulation results. 

B. RAYLEIGH 


The Rayleigh distribution is given by Eq. (19) 




2(7 


<J 


(19) 


The standard deviation in the above equation is the eommon standard deviation in 
the X and y direetions (the Rayleigh distribution requiring a eireular normal distribution 

in the x and y directions). Again, take a known (common in the x and y 

direetions) for a eireular normal dataset. For eireular normal distributions the CEP can be 
ealeulated using Eq. (14) whieh yields a CEP of 58.870 ft. Running a Monte-Carlo 
simulation to generate a eireular normal dataset for various numbers of samples and 
extracting CEP yields Table 3. 

Table 3 CEP Convergenee Data 


Number of Bombs 
Simulated 

CEP (ft) 

Run 1 

Run 2 

Run 3 

10 

62.485 

70.225 

72.979 

lo^ 

60.084 

55.831 

57.434 

10' 

58.398 

58.803 

58.136 

10" 

58.978 

59.727 

58.479 

10' 

58.918 

58.948 

58.823 

To® 

58.881 

58.883 

58.879 

10' 

58.868 

58.869 

58.865 
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c. 


SALVO FORMULA 


The salvo formula calculations are significantly more complex than the simple 
examples given for REP and CEP. However, the process remains the same. Using Eq. 
( 20 ): 


r i=i 






R. 




^ K 

2c + ^ 


V 


cr 


X J 


1-exp 


R 




r. K 

2c + ^ 

cr; 


Cl 


2cjI 


a 


X J 


( 20 ) 


With: 


n = 7 

RT=200 ft 
a^ = 150 ft 


= 50 ft 


c= 1.0 


Results in a p„(n)= 0.762. 


Again, a Monte-Carlo simulation can be performed varying the number of 
occasions to achieve convergence. The results of this simulation are shown on Table 4. 

Table 4 Salvo Eormula Convergence Data 


Number of Occasions 
Simulated 

Pn(n) 

Run 1 

Run 2 

Run 3 

10 

0.600 

0.800 

0.600 

lo^ 

0.780 

0.770 

0.720 

I(P 

0.788 

0.749 

0.787 

10^ 

0.777 

0.774 

0.778 

lo^ 

0.775 

0.772 

0.774 

10^ 

0.773 

0.773 

0.773 

10' 

0.773 

0.773 

0.773 
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It is important to remember that Eq. (20) is an analytical approximation for the 
salvo effectiveness. Using the Monte-Carlo simulation therefore can yield more accurate 
results than the analytical approximation. The approximation still gives reasonably close 
results and is only 0.011 from the converged simulation value of 0.773. However, this 
demonstrates the potential value of the Monte-Carlo simulation process by taking 
advantage of modern computing power to run enough iterations to calculate a more 
accurate result than an approximation can provide. 
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III. MAINTAINING AIMING ERROR AND BALLISTIC 
DISPERSION AS SEPARATE PARAMETERS 


In the study of statistics it is often suitable to combine multiple probabilities to 
simplify the analysis. One possible combination is taking the root-sum-square (RSS) of 
two independent standard deviations of similar distributions that are used to define the 
behavior of a total population. However, care must be taken if performing this operation 
when addressing the weapon delivery standard deviations of aiming error and ballistic 
dispersion. When attempting to calculate accuracy parameters (CEP or REP/DEP) the 
RSS simplification is appropriate. However, when calculating weapon effectiveness for a 
salvo of munitions the aiming error and ballistic dispersion must be maintained as 
separate parameters. These results are demonstrated through the use of simulation in the 
following sections. 

A. ACCURACY CALCULATIONS 

The accuracy calculations performed below are based on a simulation using an 
algorithm as outlined in Eigure 10. This algorithm assumes a circular normal distribution 
with zero mean for both error types. It is also possible to easily modify the algorithm to 
perform noncircular calculations. This is one of the significant advantages to using 
simulation practices versus analytical approaches. Many analytical approaches require 
substantial mathematical manipulation and potential approximation or numerical 
solutions to yield useful results. By creating the proper simulation routine, a very 
complex weapon accuracy model can be evaluated in a very similar fashion to the simple 
model pictured here. 
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Separate Aiming Error and Ballistic RSS of Aiming and Ballistic 


Dispersion Error Dispersion Errors 



Figure 10 Flowchart For Accuracy Simulations 

1. Single Round Scenario 

The single round accuracy results for various ratios of aiming error and ballistic 
dispersion show good correlation between the calculated values for CEP, REP, and DEP 
using both algorithms (separate errors vs. RSS error). Simulations were performed 
using 10^ occasions as outlined in Figure 10. This number of occasions allowed for 
proper convergence of the results. The results are shown below on Table 5. Results can 
also be compared to the values calculated from Eq. (12), Eq. (13) and Eq. (14). 
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Table 5 Separate vs. RSS Errors for Aeeuraey Caleulations (1 bomb per salvo) 


ERRORS 

CEP 

REP 

DEP 


Separate 

RSS 

Separate 

RSS 

Separate 

RSS 

Uaiming 50 
Obd=5 
Urss=50.3 

59.19 

59.17 

33.91 

33.90 

33.90 

33.91 

Uaiming 5 

abd=50 

Urss'^SO.S 

59.15 

59.15 

33.89 

33.88 

33.87 

33.90 

Oaiming 500 
abd=500 
nRss~707.1 

833.0 

832.6 

477.0 

476.9 

477.3 

476.8 

Claiming” 1000 
^bd^5 

Grss^I 000,01 

1,178 

1,177 

674.5 

674.2 

674.4 

675.1 

*^aiming“9 53,95 
abd=300 

Grss”1 000,01 

1,177 

1,178 

674.4 

674.4 

674.7 

674.5 


It is elear from the above table that the aiming error and ballistie dispersion error 
ean be root-sum-squared to simplify the ealeulations for the aeeuraey parameters. 
Additionally, ealeulating the RSS (or total) standard deviation from the aiming and 
ballistie dispersion standard deviations allows for the use of Eq. (12), Eq. (13) and Eq. 
(14) to easily ealeulate the aeeuraey parameters for the single round seenario. 

2. Salvo Scenario 

The salvo seenario analysis for aeeuraey ealeulations also uses the proeesses 
outlined in Eigure 10. However, the weapon loops will now be used due to the salvo 
having a greater that one number of rounds per eaeh oeeasion. Runs of 5, 10, 50 and 100 
weapons per salvo were performed with the same standard deviations as used in the 
single round seenario analysis. These results from the salvo aeeuraey parameter analysis 
ean be seen on Table 6 - Table 9. 
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Table 6 Separate vs. RSS Errors for Aeeuraey Calculations (5 bombs per salvo) 


ERRORS 

CEP 

REP 

DEP 


Separate 

RSS 

Separate 

RSS 

Separate 

RSS 

Uaiming 50 
Ubd=5 
Urss=50.3 

59.12 

59.17 

33.87 

33.89 

33.88 

33.88 

Uaiming 5 

abd=50 

Urss'^SO.S 

59.17 

59.16 

33.88 

33.92 

33.89 

33.89 

Oaiming 500 
abd=500 
nRss~707.1 

832.5 

832.8 

476.7 

477.2 

477.1 

477.0 

Claiming” 1000 
^bd^5 

Grss^I 000,01 

1,178 

1,178 

673.6 

675.1 

674.8 

674.8 

*^aiming“9 53,95 
abd=300 

Grss”1 000,01 

1,177 

1,177 

674.7 

674.2 

673.7 

673.9 


Table 7 Separate vs. RSS Errors for Accuracy Calculations (10 bombs per salvo) 


ERRORS 

CEP 

REP 

DEP 


Separate 

RSS 

Separate 

RSS 

Separate 

RSS 

*^aiming“50 

*^bd“5 

C^RSS“50,3 

59.11 

59.18 

33.86 

33.90 

33.90 

33.91 

Uaiming 5 
abd=50 
Urss=50.3 

59.19 

59.15 

33.9 

33.89 

33.90 

33.90 

Uaiming 500 
abd=500 
Urss“707.1 

832.7 

832.6 

477.1 

477.0 

477.1 

477.1 

*^aiming~ 1000 
*^bd“5 

Grss“1 000,01 

1,177 

1,178 

673.3 

674.5 

675.3 

674.5 

*^aimmg~953,95 
Gbd=300 
Grss~1 000,01 

1,177 

1,178 

674.5 

674.5 

674.2 

674.7 
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Table 8 Separate vs. RSS Errors for Aecuracy Calculations (50 bombs per salvo) 


ERRORS 

CEP 

REP 

DEP 


Separate 

RSS 

Separate 

RSS 

Separate 

RSS 

Uaiming 50 

nbd=5 

O'rss~50.3 

59.15 

59.16 

33.86 

33.89 

33.84 

33.89 

Uaiming 5 

nbd=50 

f^RSS~50.3 

59.16 

59.16 

33.89 

33.89 

33.89 

33.89 

Uaiming 500 
abd=500 
f^RSS“207.1 

832.8 

832.4 

477.4 

476.9 

476.9 

476.9 

*^aimmg” 1000 

Grss^IOOO.01 

1,179 

1,177 

674.2 

674.4 

676.3 

674.4 

*^aiming~9 53,95 
abd=300 
Grss^IOOO.OI 

1,177 

1,178 

675.3 

674.5 

674.1 

674.6 


Table 9 Separate vs. RSS Errors for Accuracy Calculations (100 bombs per salvo) 


ERRORS 

CEP 

REP 

DEP 


Separate 

RSS 

Separate 

RSS 

Separate 

RSS 

Uaiming 50 

nbd=5 

nRss=50.3 

59.20 

59.17 

33.91 

33.89 

33.91 

33.89 

Uaiming 5 
abd=50 
nRss'^50.3 

59.17 

59.16 

33.90 

33.90 

33.90 

33.89 

Oaiming 500 
abd=500 
Urss~707.1 

832.6 

832.5 

476.9 

476.9 

477.1 

476.9 

*^aiming” 1000 

Grss^I 000.01 

1,177 

1,177 

673.7 

674.5 

673.6 

674.5 

*^aiming“953.95 

abd=300 

Grss”1 000.01 

1,178 

1,177 

674.8 

674.5 

674.7 

674.5 
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It is interesting to note that on Table 5 thru Table 9 the two methods for 
ealeulating the aeeuraey parameters both show similar results for the various seenarios 
regardless of the number of weapons per salvo. This demonstrates that for aeeuraey 
parameter ealeulation it is appropriate to RSS the aiming error and ballistie dispersion 
error standard deviations for both single round and salvo seenarios. 

B. WEAPON EFFECTIVENESS CALCULATIONS 

For weapon effectiveness simulations the procedure outlined in Figure 11 was 
used. This algorithm can be easily modified to allow for a more complex weapon/target 
lethal area than is available with analytical approaches. 
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Separate Aiming Error and Ballistic 
Dispersion Error 


RSS of Aiming and Ballistic 
Dispersion Errors 




Figure 11 Flow Chart for Weapon Effectiveness Simulations 

1. Single Round Scenario 

The single round scenario will provide reasonable results using both of the 
methods outlined in Figure 11. The Mean Area of Effectiveness was set to 7854 
corresponding to a lethal radius of 50 ft. Simulation results can be seen in Table 10. 
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An analytical approach can also be used to eheek the single round seenario 
algorithms. Using the CDF for the Rayleigh distribution will yield the probability that a 
given round will fall within the radius R. If R is set to the lethal radius and the standard 
deviation is set to the RSS value of the given aiming and ballistie dispersion errors the 
CDF value will equal the probability that a given round will land within the lethal radius 
of the weapon. This is preeisely the weapon effeetiveness parameter of interest for the 
single round seenario. Knowing that the Rayleigh CDF is: 




( 21 ) 


Where: R = Weapon Lethal Radius = 50 ft 


cr = 


2 2 
*^aiming ' ^bd 


Table 10 Separate vs. RSS Errors for Weapon Effeetiveness (1 bomb per salvo) 


ERRORS 

Pn(: 

1) 

Eq (21) 


Separate 

RSS 

Analytical Sol. 

Uaiming 50 
Ubd=5 
Urss'^SO.S 

0.39 

0.39 

0.39 

Uaiming 5 

abd=50 

Urss="50.3 

0.39 

0.39 

0.39 

Uaiming 500 
abd=500 
Urss“707.1 

0.0025 

0.0025 

0.0025 

*^aiming“ 1000 
^bd^5 

Grss”1000,01 

0.0012 

0.0012 

0.0012 

aiming“9 53.95 

‘^bd^300 

Grss”1000,01 

0.0012 

0.0012 

0.0012 


It is elear from the results in Table 10 that either method is suitable for 
determining the weapon system effeetiveness for a single round seenario. 
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2 . 


Salvo Scenario 


The salvo scenario analysis for accuracy calculations uses the processes outlined 
in Figure 11. However, the weapon loops will now be used due to the salvo having a 
greater that one number of rounds per occasion. Runs of 5, 10, 50 and 100 weapons per 
salvo were performed with the same standard deviations as used in the single round 
scenario analysis. These results from the salvo weapon effectiveness analysis can be seen 
on Table 11. 


Table 11 Separate vs. RSS Errors for Weapon Effectiveness (5,10,50,100 bombs per salvo) 


ERRORS 

Pn(5) 

Pn(lO) 

Pn(50) 

Pn(lOO) 


Separate 

RSS 

Separate 

RSS 

Separate 

RSS 

Separate 

RSS 

Uaiming 50 

nbd=5 

f^RSS~50.3 

0.46 

0.92 

0.48 

1.00 

0.52 

1.00 

0.54 

1.00 

Uaiming 5 

abd=50 

Urss^SO.S 

0.92 

0.92 

0.99 

0.99 

1.00 

1.00 

1.00 

1.00 

Uaiming 500 
abd=500 
O'RSS^^Oy.l 

0.012 

0.012 

0.025 

0.025 

0.12 

0.12 

0.21 

0.23 

*^aiming~ 1000 
0^bd^5 

(Jrss~ 1000.01 

0.0016 

0.0062 

0.0016 

0.012 

0.0021 

0.061 

0.0018 

0.12 

aiming"953,95 

abd=300 

<Jrss" 1000.01 

0.0062 

0.0062 

0.012 

0.012 

0.053 

0.059 

0.092 

0.11 


Table 11 demonstrates the fundamental reason that the aiming error and ballistic 
dispersion error standard deviations must be treated separately for the purposes of 
weapon effectiveness calculations. The individual weapons of a given salvo are all 
influenced by the same aiming error. This results in the weapon impact location for each 
weapon in the salvo being dependent on a constant aiming error. Eor a given salvo to be 
considered successful in killing the target at least one bomb of that salvo must impact 
within the lethal area. The single round separate error scenario however, will count each 
weapon inside the lethal area as a kill and consequently result in a similar effectiveness 
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parameter as that of the RSS weapon effeetiveness. To eheek this assumption, an 
algorithm ean be setup to eount eaeh bomb that lands inside the lethal area during the 
separate algorithm operation. Count all the bombs that fall inside the lethal area and 
divide by the total number of bombs dropped to yield the ehanee of a given bomb landing 
inside the lethal area. Then, this number must be powered up to provide the ineorreet 
total salvo weapon effeetiveness that will eorrespond to the ineorreet value given by the 
RSS algorithm. This fundamental differenee results in the RSS algorithm providing 
ineorreetly optimistie weapon effeetiveness for weapon salvos. For the seenarios where 
the aiming error is signifieantly large than the ballistie dispersion error the separate vs. 
RSS effeetiveness values are signifieantly different. This error is a funetion of Oaiming, 
Obd, and lethal radius. Clearly, for proper ealeulation of weapon effeetiveness, the aiming 
error and ballistie dispersion errors must be kept separate and the proeedure on the left 
side of Figure 11 should be used. 
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IV. SALVO EFFECTIVENESS CALCULATIONS: SIMULATION 
VS. ANALYTICAL APPROXIMATIONS 


It was observed in Chapter II that the simulation for salvo weapon effeetiveness 
eould yield more aeeurate results than the analytieal approximations historieally used. 
Chapter IV will provide additional detail regarding these approximation defieieneies. 

A. CIRCULAR TARGET 

The following eharts eompare simulation runs with Eq. (16) and Eq. (17) salvo 
formula approximations. Two sets of errors were evaluated. One set with aiming error 
of 50 and ballistie dispersion of 5 (as investigated previously). The other set inereased 
the ballistie dispersion and uses aiming error of 50 and ballistie dispersion of 25. 
Weapon effeetiveness for both of these error sets was ealeulated for multiple weapon 
lethal areas (5, 10, 20, 30, 40, 50, 60). The results are shown below on Table 12 thru 
Table 25. 


Table 12 Salvo Sim. vs. Approximation: (aaiming=50; abd=5; Eethal Radius 5) 


Lethal Radius 

5 


#/Salvo 

Sim Circle 

Eq. 16 

Eq. 17 

Delta: Sim-Eq 16 

Delta: Sim-Eq 17 

1 

0.005 

0.005 

0.005 

0.000 

0.000 

10 

0.025 

0.027 

0.027 

-0.002 

-0.002 

20 

0.033 

0.037 

0.037 

-0.003 

-0.003 

30 

0.038 

0.042 

0.042 

-0.004 

-0.004 

40 

0.041 

0.046 

0.046 

-0.005 

-0.005 

50 

0.044 

0.049 

0.050 

-0.006 

-0.006 

60 

0.046 

0.052 

0.052 

-0.006 

-0.006 

70 

0.048 

0.054 

0.054 

-0.006 

-0.006 

80 

0.049 

0.056 

0.056 

-0.008 

-0.008 

90 

0.050 

0.058 

0.058 

-0.007 

-0.008 

100 

0.052 

0.059 

0.059 

-0.007 

-0.007 
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Table 13 Salvo Sim. vs. Approximation: (aaiming=50; abd=5; Lethal Radius 10) 


Lethal Radius 

10 


#/Salvo 

Sim Circle 

Eq 16 

Eq 17 

Delta: Sim-Eq 16 

Delta: Sim-Eq 17 

1 

0.020 


0.019 

0.001 

0.000 

10 

0.055 

0.072 

0.072 

-0.016 

-0.017 

20 

0.066 

0.090 

0.091 

-0.024 

-0.024 

30 

0.073 

0.100 

0.101 

-0.028 

-0.029 

40 

0.078 

0.108 

0.109 

-0.030 

-0.031 

50 

0.080 

0.114 

0.115 

-0.034 

-0.035 

60 

0.082 

0.119 

0.120 

-0.037 

-0.038 

70 

0.085 

0.120 

0.124 

-0.034 

-0.038 

80 

0.087 

-3.42E-01 

0.123 

0.430 

-0.035 

90 

0.090 

-5.73E+01 

-6.13E-01 

5.74E+01 

7.04E-01 

100 

0.090 

-6.07E+03 

-1.33E+02 

6.07E+03 

1.33E+02 


Table 14 Salvo Sim. vs. Approximation: (aaiming=50; abd=5; Lethal Radius 20) 


Lethal Radius 

20 


#/Salvo 

Sim Circle 

Eq 16 

Eq 17 

Delta: Sim-Eq 16 

Delta: Sim-Eq 17 

1 

0.076 

0.071 

0.073 

0.005 

0.003 

10 

0.138 

0.211 

0.219 

-0.074 

-0.081 

20 

0.153 

0.256 

0.264 

-0.103 

-0.112 

30 

0.163 

0.281 

0.290 

-0.118 

-0.127 

40 

0.169 

0.298 

0.308 

-0.130 

-0.139 

50 

0.173 

0.311 

0.322 

-0.139 

-0.149 

60 

0.176 

0.298 

0.330 

-0.122 

-0.154 

70 

0.179 

-1.36E+01 

-2.46E+00 

1.38E+01 

2.64E+00 

80 

0.181 

-3.08E+03 

-6.13E+02 

3.08E+03 

6.14E+02 

90 

0.182 

-5.17E+05 

-8.60E+05 

5.17E+05 

8.60E+05 

100 

0.186 

5.44E+07 

-2.31E+08 

-5.44E+07 

2.31E+08 
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Table 15 Salvo Sim. vs. Approximation: (aaiming=50; abd=5; Lethal Radius 30) 


Lethal Radius 

30 


#/Salvo 

Sim Circle 

Eq 16 

Eq 17 

Delta: Sim-Eq 16 

Delta: Sim-Eq 17 

1 

0.163 

0.140 

0.151 

0.022 

0.011 

10 

0.244 

0.380 

0.406 

-0.136 

-0.162 

20 

0.262 

0.448 

0.476 

-0.186 

-0.214 

30 

0.273 

0.485 

0.514 

-0.212 

-0.241 

40 

0.281 

0.510 

0.540 

-0.229 

-0.259 

50 

0.285 

0.528 

0.559 

-0.243 

-0.274 

60 

0.288 

0.592 

0.592 

-0.304 

-0.304 

70 

0.293 

-5.44E+01 

-6.82E+01 

5.47E+01 

6.85E+01 

80 

0.296 

1.51E+04 

-8.73E+03 

-1.51E+04 

8.73E+03 

90 

0.298 

-4.75E+07 

-1.47E+07 

4.75E+07 

1.47E+07 

100 

0.301 

-2.88E+10 

-l.OlE+10 

2.88E+10 

l.OlE+10 


Table 16 Salvo Sim. vs. Approximation: (aaiming=50; abd=5; Lethal Radius 40) 


Lethal Radius 

40 


#/Salvo 

Sim Circle 

Eq 16 

Eq 17 

Delta: Sim-Eq 16 

Delta: Sim-Eq 17 

1 


■BM 

0.241 

0.057 

0.030 

10 



0.587 


-0.223 

20 


MEm 

0.668 


-0.286 

30 

0.395 

0.656 

0.708 

-0.261 

-0.314 

40 

0.401 

0.683 

0.734 

-0.282 

-0.334 

50 

0.406 

0.702 

0.753 

-0.296 

-0.347 

60 

0.408 

0.797 

0.930 

-0.389 

-0.522 

70 

0.412 

5.18E+01 

-l.OlE+02 

-5.14E+01 

1.02E+02 

80 

0.416 

1.46E+05 

-5.11E+04 

-1.46E+05 

5.11E+04 

90 

0.419 

-1.18E+08 

-6.46E+07 

1.18E+08 

6.46E+07 

100 

0.420 

-1.81E+11 

-6.06E+10 

1.81E+11 

6.06E+10 
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Table 17 Salvo Sim. vs. Approximation: (aaiming=50; abd=5; Lethal Radius 50) 


Lethal Radius 

50 


#/Salvo 

Sim Circle 

Eq 16 

Eq 17 

Delta: Sim-Eq 16 

Delta: Sim-Eq 17 

1 

0.390 

0.282 

0.331 

0.108 

0.059 

10 

0.485 

0.658 

0.733 

-0.173 

-0.248 

20 

0.503 

0.739 

0.809 

-0.236 

-0.307 

30 

0.513 

0.778 

0.844 

-0.265 

-0.331 

40 

0.519 

0.802 

0.865 

-0.283 

-0.346 

50 

0.524 

0.819 

0.880 

-0.295 

-0.356 

60 

0.529 

0.790 

1.19E+00 

-0.261 

-0.665 

70 

0.535 

-1.87E+02 

-3.79E+02 

1.88E+02 

3.80E+02 

80 

0.535 

9.30E+04 

1.51E+05 

-9.30E+04 

-1.51E+05 

90 

0.538 

-3.42E+08 

-1.76E+08 

3.42E+08 

1.76E+08 

100 

0.540 

-1.72E+11 

-8.75E+10 

1.72E+11 

8.75E+10 


Table 18 Salvo Sim. vs. Approximation: (aaiming=50; abd=5; Lethal Radius 60) 


Lethal Radius 

60 


#/Salvo 

Sim Circle 

Eq 16 

Eq 17 

Delta: Sim-Eq 16 

Delta: Sim-Eq 17 

1 

0.510 

0.340 

0.416 

0.169 

0.093 

10 

0.597 

0.748 

0.838 

-0.151 

-0.241 

20 

0.613 

0.824 

0.899 

-0.211 

-0.287 

30 

0.622 

0.858 

0.924 

-0.236 

-0.302 

40 

0.630 

0.878 

0.938 

-0.249 

-0.308 

50 

0.636 

0.892 

0.949 

-0.255 

-0.313 

60 

0.639 

1.75E+00 

1.91E+00 

-l.llE+00 

-1.27E+00 

70 

0.642 

4.83E+01 

-1.91E+02 

-4.77E+01 

1.91E+02 

80 

0.645 

4.77E+05 

1.03E+04 

-4.77E+05 

-1.03E+04 

90 

0.647 

-1.37E+08 

-3.85E+08 

1.37E+08 

3.85E+08 

100 

0.648 

9.58E+10 

-2.28E+11 

-9.58E+10 

2.28E+11 
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Table 19 Salvo Sim. vs. Approximation: (aaiming=50; abd=25; Lethal Radius 5) 


Lethal Radius 

5 


#/Salvo 

Sim Circle 

Eq 16 

Eq 17 

Delta: Sim-Eq 16 

Delta: Sim-Eq 17 

1 

0.004 

0.004 

0.004 

0.000 

0.000 

10 

0.038 

0.038 

0.038 

0.000 

0.000 

20 

0.073 

0.072 

0.072 

0.001 

0.001 

30 

0.104 

0.103 

0.103 

0.001 

0.001 

40 

0.131 

0.130 

0.131 

0.001 

0.001 

50 

0.156 

0.155 

0.156 

0.001 

0.001 

60 

0.179 

0.178 

0.179 

0.001 

0.001 

70 

0.198 

0.199 


-0.001 

-0.002 

80 

0.219 

0.218 

wm 

0.001 

0.001 

90 

0.238 

0.236 


0.002 

0.001 

100 

0.252 

0.252 

0.252 

0.000 

0.000 


Table 20 Salvo Sim. vs. Approximation: (aaiming=50; abd=25; Lethal Radius 10) 


Lethal Radius 

10 


#/Salvo 

Sim Circle 

Eq 16 

Eq 17 

Delta: Sim-Eq 16 

Delta: Sim-Eq 17 

1 

0.016 

0.016 

0.016 


0.000 

10 

0.131 

0.131 

0.132 


0.000 

20 

0.220 

0.220 

0.222 

0.000 

-0.001 

30 

0.283 

0.283 

0.285 

0.000 

-0.002 

40 

0.331 

0.331 

0.333 

0.001 

-0.002 

50 

0.368 

0.367 

0.370 

0.001 

-0.002 

60 

0.391 

0.396 

0.399 

-0.005 

-0.008 

70 

0.418 

0.420 

0.423 

-0.002 

-0.005 

80 

0.437 

0.440 

0.443 

-0.003 

-0.006 

90 

0.453 

0.457 

0.460 

-0.004 

-0.007 

100 

0.468 

0.472 

0.475 

-0.004 

-0.007 
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Table 21 Salvo Sim. vs. Approximation; (aaiming=50; abd=25; Lethal Radius 20) 


Lethal Radius 

20 


#/Salvo 

Sim Circle 

Eq 16 

Eq 17 

Delta: Sim-Eq 16 

Delta: Sim-Eq 17 

1 

0.061 

0.058 

0.060 

0.003 

0.001 

10 

0.345 

0.343 

0.352 

0.002 

-0.008 

20 

0.459 

0.463 

0.475 

-0.004 

-0.016 

30 

0.515 

0.526 

0.539 

-0.012 

-0.025 

40 

0.555 

0.567 

0.580 

-0.012 

-0.025 

50 

0.581 

0.596 

0.610 

-0.015 

-0.029 

60 

0.601 

0.619 

0.632 

-0.017 

-0.031 

70 

0.620 

0.637 



-0.030 

80 

0.635 

0.652 



-0.030 

90 

0.646 

0.664 


-0.018 

-0.032 

100 

0.656 

0.675 

0.689 

-0.019 

-0.033 


Table 22 Salvo Sim. vs. Approximation; (aaiming=50; abd=25; Lethal Radius 30) 


Lethal Radius 

30 


#/Salvo 

Sim Circle 

Eq 16 

Eq 17 

Delta: Sim-Eq 16 

Delta: Sim-Eq 17 

1 

0.135 

0.118 

0.126 

0.017 

0.009 

10 

0.507 

0.509 

0.536 

-0.001 

-0.028 

20 

0.605 

0.621 

0.650 

-0.016 

-0.045 

30 

0.652 

0.676 

0.705 

-0.024 

-0.053 

40 

0.686 

0.710 

0.738 

-0.024 

-0.053 

50 

0.707 

0.734 

0.762 

-0.027 

-0.055 

60 

0.723 

0.752 

0.780 

-0.029 

-0.056 

70 

0.735 

0.767 

0.794 

-0.031 

-0.058 

80 

0.747 

0.778 

0.805 

-0.031 

-0.058 

90 

0.755 

0.788 

0.815 

-0.033 

-0.059 

100 

0.765 

0.796 

0.822 

-0.032 

-0.057 
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Table 23 Salvo Sim. vs. Approximation; (aaiming=50; abd=25; Lethal Radius 40) 


Lethal Radius 

40 


#/Salvo 

Sim Circle 

Eq 16 

Eq 17 

Delta: Sim-Eq 16 

Delta: Sim-Eq 17 

1 

0.226 

0.184 

0.204 

0.041 

0.022 

10 

0.630 

0.633 

0.681 

-0.003 

-0.051 

20 

0.715 

0.734 

0.781 

-0.019 

-0.065 

30 

0.756 

0.781 

0.825 

-0.025 

-0.069 

40 

0.781 

0.809 

0.851 

-0.029 

-0.070 

50 

0.797 

0.829 

0.868 

-0.032 

-0.071 

60 

0.810 

0.844 

0.881 

-0.034 

-0.071 

70 

0.819 

0.855 

0.891 

-0.036 

-0.072 

80 

0.828 

0.865 

0.899 

-0.037 

-0.071 

90 

0.835 

0.877 

0.861 

-0.041 

-0.026 

100 

0.842 

0.933 

-5.13E+00 

-0.092 

5.97E+00 


Table 24 Salvo Sim. vs. Approximation; (aaiming=50; abd=25; Lethal Radius 50) 


Lethal Radius 

50 


#/Salvo 

Sim Circle 

Eq 16 

Eq 17 

Delta: Sim-Eq 16 

Delta: Sim-Eq 17 

1 

0.331 

0.249 

0.286 

0.082 

0.045 

10 

0.732 

0.725 

0.792 

0.007 

-0.060 

20 

0.801 

0.815 

0.872 

-0.015 

-0.072 

30 

0.832 

0.854 

0.905 

-0.023 

-0.073 

40 

0.849 

0.877 

0.923 

-0.028 

-0.074 

50 

0.864 

0.893 

0.935 

-0.029 

-0.071 

60 

0.872 

0.904 

0.943 

-0.032 

-0.071 

70 

0.881 

0.909 

0.943 

-0.028 

-0.063 

80 

0.886 

0.913 

0.565 

-0.026 

0.321 

90 

0.893 

-2.80E+01 

1.79E+01 

2.89E+01 

-1.70E+01 

100 

0.896 

3.83E+03 

1.28E+03 

-3.83E+03 

-1.28E+03 
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Table 25 Salvo Sim. vs. Approximation; (aaiming=50; abd=25; Lethal Radius 60) 


Lethal Radius 

60 


#/Salvo 

Sim Circle 

Eq 16 

Eq 17 

Delta: Sim-Eq 16 

Delta: Sim-Eq 17 

1 

0.438 

0.306 

0.365 

0.132 

0.072 

10 

0.812 

0.793 

0.871 

0.019 

-0.059 

20 

0.865 

0.871 

0.931 

-0.006 

-0.066 

30 

0.888 

0.903 

0.953 

-0.016 

-0.065 

40 

0.902 

0.921 

0.964 

-0.019 

-0.062 

50 

0.911 

0.933 

0.971 

-0.022 

-0.060 

60 

0.919 

0.942 

0.975 

-0.023 

-0.056 

70 

0.925 

0.824 

0.853 

0.100 

0.072 

80 

0.929 

0.457 

-4.68E+00 

0.473 

5.610 

90 

0.932 

-3.96E+03 

-4.90E+03 

3.96E+03 

4.90E+03 

100 

0.934 

1.26E+05 

-5.47E+05 

-1.26E+05 

5.47E+05 


Two speeific results can be observed from the tables above. First, it is clear that 
the approximations in Eq. (16) and Eq. (17) will break down at high values of number of 
bombs per salvo. This behavior is caused by the “battle of big alternating binomial 
coefficients” that results for high values of n [2]. The sudden breakdown of weapon 
effectiveness values calculated from Eq. (16) and Eq. (17) demonstrates one of the risks 
of using such approximations. 

Second, the approximations performed rather poorly for Table 12 thru Table 18 
when the ballistic dispersion was significantly smaller than aiming error. This 
demonstrates a more significant shortcoming of the approximations. This scenario is 
becoming more and more common as guided weapons become extremely accurate and 
self-designation by the fighter/bomber can cause a relatively large aiming error. 

Einally, with the capability of modern computing systems and a significant level 
of uncertainty regarding the outputs of Eq. (16) and Eq. (17), it is recommended to use 
the separate error salvo weapon effectiveness procedure as detailed in Eigure 11 
whenever possible. 
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B. SQUARE TARGET 


The data in the previous seetion used a eireular MAE. It is possible that a weapon 
effeetive area is more reetangular in nature. More of the weapon effeetiveness is foeused 
out of the side of the munition easing. If a bomb impaets with a shallow angle the area of 
effeetiveness is longer in the defleetion direetion and shorter in the range direetion. In 
the absenee of a salvo weapon effeetiveness equation/approximation the previous weapon 
effeetiveness proeedure from Figure 11 ean be altered to eheek if the weapon impaets a 
reetangular area eentered on the target. This again demonstrates the versatility of using 
the simulation proeedure to ealeulate weapon effeetiveness parameters. 

The following results show the weapon effeetiveness values for a eirele of radius 
25 eompared to various aspeet ratio reetangles of the same area. Length to width ratios 
of 1, 0.5, 0.25, 0.125, and 0.0625 will be evaluated for the aaiming=50; abd=25; Lethal 
Radius 20 seenario. 

Table 26 Cireular vs. Reetangle Weapon Effeetiveness Area Salvo Simulation 


Oaimine=50; abd=25; Lethal Rad 

lius 20 


Ratio: 

1 

0.5 

0.25 

0.125 

0.0625 

#/Salvo 

Circle 

Square 

Rectangle 

Rectangle 

Rectangle 

Rectangle 

1 

0.062 

0.062 

0.062 

0.059 

0.056 

0.050 

10 

0.346 

0.509 

0.504 

0.488 

0.452 

0.397 

20 

0.458 

0.616 

0.617 

0.610 

0.596 

0.553 

30 

0.517 

0.667 

0.670 

0.669 

0.660 

0.628 

40 

0.554 

0.700 

0.701 

0.704 

0.694 

0.671 

50 

0.581 

0.722 

0.723 

0.725 

0.718 

0.701 

60 

0.602 

0.739 

0.743 

0.743 

0.737 

0.722 

70 

0.621 

0.755 

0.755 

0.755 

0.752 

0.737 

80 

0.633 

0.763 

0.765 

0.769 

0.765 

0.751 

90 

0.646 

0.774 

0.775 

0.778 

0.775 

0.760 

100 

0.658 

0.782 

0.782 

0.784 

0.782 

0.771 
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The weapon effeetiveness comparisons shown in Table 26 detail a scenario for a 
circular distribution of impacts with a rectangular MAE for the weapon. The rectangular 
MAE results in slightly higher weapon effectiveness than the circular MAE of equivalent 
area. Also of interest are the slightly smaller effectiveness parameters that result from the 
reduced MAE aspect. This implies, for this scenario, as the rectangular MEA becomes 
thinner and longer the weapon effectiveness is slightly reduced. With minor 
modifications these techniques can also be used to calculate weapon effectiveness for 
non-circular impact distributions with any weapon/target MAE shape desired. 
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V. ATTEMPTS TO EXTRACT ERROR TYPES EROM IMPACT 

DATA 


Often during weapon delivery performanee testing the primary data gathered are 
the range and defleetion miss distanees. Having the weapon impaet miss distanees and 
an estimation of the type of expeeted distribution one would hope to be able to extract the 
aiming error and ballistic dispersion error values that created the impact locations. 
Unfortunately, for a normal distribution this extraction is ambiguous due to the RSS 
value of the aiming and ballistic dispersion errors allowing for an infinite number of 
combinations for the same distribution. However, there is still value added in attempting 
to properly characterize the weapon impact distribution for the purpose of weapon 
effectiveness calculations. Complicating this process is the fact that guided weapon 
distributions are typically not normal distributions. Therefore, a method to approximate 
an appropriate distribution for a non-normal dataset must be introduced. 

A. ALGORITHM DESCRIPTION 

This process starts by assuming a particular form for the weapon miss distance 
CDF. For the purposes of this discussion, it will be assumed that the non-normal 
distribution is an expression derived from linear combinations of other known 
distributions weighted accordingly. These distributions will be discussed in more detail 
in Section B. Once a distribution is defined, the impact data needs to be rank ordered and 
a CDF needs to be calculated from the impact data. Once this empirically derived CDF is 
determined a least-squares-fit to the assumed CDF is performed using MATLAB. The 
parameters being used to accomplish the curve fit are the standard deviations of the CDF 
and the weighting factors. Once the weighting factors and standard deviations are 
known, the weapon effectiveness calculations can be performed. See Figure 12 for a 
diagram of this process. See Section B for examples. [4] 
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Figure 12 Flow Chart for Extracting Weapon Effectiveness Erom Impact Data 

B. SINGLE ROUND SCENARIO 


I. Double Normal Approximation to Non-Normal Dataset 


Eor an assumed double normal dataset the distribution would be defined as 
follows. Knowing that the univariate normal PDE is: 


fix) = 



2a^ 


And the univariate normal CDE is: 


( 22 ) 
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( 23 ) 


F(x)= f— ^ dx = — 1 + erf 
Where: z = ——— 

<j 

Creating the double normal CDF in the range direetion yields: 

CDF™. l + 

1 r ( 1 r ( V 

Where: w^j + w ^2 = 1 



(24) 



As a test of this algorithm, a dataset will be ereated using Eq. (25). Setting 
w^j =0.7, cr_^j =30, cr ^2 =5 over the range [-100:2:100] and adding a small error souree 
to randomize the dataset results in the dataset shown in Figure 13. 
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Figure 13 Artificially Created Double Normal CDF ( = 0.3, = 30, = 5 ) 

The next step is to curve fit the data to extract vv_^i, cr_^i, crj. 2 . Using initial 
conditions of: 

w,i=0.1 =0.317 

cr_^j = 5 results in curve fitting estimates of: cr^j = 5.693 
= 40 = 30.407 

The curve created by these estimates can be seen in comparison with the test 
dataset on Figure 14. 









B3 - eoo 

fO - 





o' 









o 

^ O 









o 

0 








c 

o 








o 









^ “O 

o 

o 









o 

o 

o 

_ 








d 

cP 








yC>S§St3X 









44 

























Miss distance (ft) 


Figure 14 Comparison of CDFs; Double Normal Dataset vs. Curve Fitting 

It is elear for this test ease the algorithm eurve fits the dataset reasonably well. It 
is interesting to note that cr_^j and cr _^2 flipped due to the procedure converging to a 

of 0.317 instead of the actual value of 0.7. This potential swapping of values between 

the two standard deviations will not have an effect of future weapon effectiveness 
calculations as the weights associated with them are also swapped. A similar process can 
be performed for the deflection direction. 
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a. Double Normal Dataset Weapon Effectiveness Calculations 

Once cr^j, been extracted from the non-normal 

datasets in the range and deflection direction the weapon effectiveness calculations can 
be performed. Again, a MAE will be required to perform the weapon effectiveness 
simulations. This MAE will be assumed to be rectangular with length (in the range 
direction) Eet and width (in the deflection direction) Wet- 

The range direction calculations for weapon effectiveness can again be 
performed using Monte-Carlo simulations. Each range standard deviation (cr^j and cr^j) 
can be simulated individually to determine the weapon effectiveness resulting from each 
parameter. These Monte-Carlo simulations will be performed as outlined on the right 
side of Eigure 11. To determine if an individual weapon is a kill the impact location is 
compared to Eet/2. If inside Eet/2 the weapon is considered a kill. These two 
simulations result what is called a single sortie probability of damage (SSPD). In the 
range direction the two SSPDs will be combined using Eq. (27) to yield a total SSPD in 
the range direction. 

+ (1 - * SSPD^2 (27) 

hollowing a similar procedure for the deflection direction using WET/2 to 
determine a kill results in Eq. (28). 

* SSPD^, + (1 - w^,) * SSPD^2 (28) 

To complete the weapon effectiveness calculation the total SSPD is 
defined as Eq. (29). 

SSPD^o^.,=SSPD^ *SSPD^ (29) 

As an example, using the previous curve fitted values of 
=0.317, cr_^j =5.693, and =30.407 and running the Monte Carlo simulations for 
the range direction assuming Eet=20 and Wet=20 results in; 
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SSPD^^ = 0.922 
SSPD^^ = 0.257 

= ^.r'^SSPD^.+il-wJ^SSPD^, 

SSPD^ =0.317 *0.922+ (1-0.317) *0.257 = 0.47 

^Total ^ ^ 


Assuming the distribution is circular results in: 


SSPD. 


TOTAL 


= SSPD^ 


-SSPD 


yTOTAL 


= SSPD^ 


SSPD^^^^^ = 0.47" = 0.22 


(31) 


2. Double Rayleigh Approximation to Non-Normal Dataset 


The procedure is similar to the one describe for a double normal distribution but 
instead a linear combination of Rayleigh distributions is used. Knowing the Rayleigh 
PDF is: 


f(r) = ^P 


a 


And the Rayleigh CDF is: 


2a‘^ 




(32) 


(33) 


F(R) = l-eL 

Substituting this expression into the correct terms of Eq. (25) yields the double 
Rayleigh distribution in Eq. (34). 


(34) 


As a test of this algorithm, a dataset will be created using Eq. (34). Setting 
=0.7, cr_^j =30, cr ^2 =5 over the range [0:2:100] and adding a small error source to 
randomize the dataset results in the dataset shown in Eigure 15. 
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Figure 15 Artificially Created Double Rayleigh CDF (= 0.3, cr^j = 30, = 5 ) 

The next step is to curve fit the data to extract vv^i,cr^i,0'r2 • Using initial 
conditions of; 

w^j=0.1 w^j= 0.300 

cr^j = 5 results in curve fitting estimates of cr^j = 4.794 
a^,=40 =29.157 

The curve created by these estimates can be seen in comparison with the test 
dataset on Figure 16. 
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Figure 16 Comparison of CDFs; Double Rayleigh Dataset vs. Curve Fitting 

Again, similarly to the double normal distribution extraction tests, the double 
Rayleigh procedure produced a reasonable curve fit and extracted values for the standard 
deviation and weighting value close to the original values. These values can now be used 
to perform weapon effectiveness calculations. 

a. Double Rayleigh Dataset Weapon Effectiveness Calculations 

The weapon effectiveness calculations for the double Rayleigh distribution 
are more straightforward than those for the double normal distribution because there is 
only one axis to analyze. Therefore, a simulation is made for each standard deviation to 
generate a radial miss distance and compare this with a predefined lethal radius. Each of 
these simulations will output an SSPDr. Once both SSPDri and SSPDr2 are simulated the 
total SSPD can be calculated as shown in Eq. (35). 
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SSPD^ = * SSPD^^ + (1 - ) * SSPD^^ 

^Totai n V r i / rz 


( 35 ) 


As an example, using the previous eurve fitted values of 
= 0.300, cr^j = 4.794, and a ^2 ~ 29.157 and running the Monte Carlo simulations for 
the radial direetion assuming Weapon Lethal Radius=l 1.28 results in; 

SSPD^^ = 0.937 
SSPD^2 = 0023 

SSPD,^,2^, = w^,*SSPD^,+{\-w^,rSSPD^2 
SSPD^o^^2. = 0.300 * 0.937 + (1 - 0.300) * 0.073 = 0.33 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


From the numerical investigations presented above, the following can be 
concluded; 

1. Extraction of aiming error and ballistic dispersion error from single 
weapon impact test data is not possible due to an infinite number of 
combinations of errors that will yield the same distribution. 

2. If salvo weapon impact data is available, aiming error and ballistic 
dispersion error extraction would be possible. 

3. Aiming error and ballistic dispersion error can be combined using a root- 
sum-square for the purposes of accuracy parameter calculations and 
simulations. 

4. The aiming error and ballistic dispersion cannot be combined using a root- 
sum-square for the purposes of weapon effectiveness calculations. 

5. Monte-Carlo simulation for salvo weapon effectiveness can provide more 
accurate results than salvo equation approximations 

6. For single round scenarios it is possible to extract two standard deviations 
with associated weighting values and calculate weapon effectiveness for a 
single round. 

The following are recommendations for future research; 

1. For salvo weapon impact data error extraction, investigate the required 
number of bombs per salvo and number of salvos required to yield 
accurate aiming error and ballistic dispersion error. 

2. Investigate the use of MATFAB Symbolic Toolbox to provide numerical 
answers to derived analytical expressions for weapon effectiveness. 
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3. Investigate the potential of MATLAB Statistics Toolbox to identify 
dataset characteristics to create more accurate non-normal distribution 
functions to be used for standard deviation extraction. 

4. Continue researching the sensitivity of the salvo weapon effectiveness 
result by varying the ratio of aiming error to ballistic dispersion error, the 
number of weapons per salvo, and the lethal radius. 
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APPENDIX 


MATLAB CODE 


A. CHAPTER II CODE 

I. Code for Table 2 


clear 

sigmax=50 

%REP calculation 
bombnumber=l0 00000; 
x=ones(bombnumber,10); 
for k=l:1:10 
1 = 1 ; 

while i<=bombnumber 
X(1,k)=randn*sigmax; 
1 = 1 + 1 ; 
end 
end 

REP=median(abs(x)) 


2. Code for Table 3 


clear 
sigmax=50; 
sigmay=50; 

%REP calculation 
bombnumber=l0 000000; 
x=ones(bombnumber,1); 
y=ones(bombnumber,1); 
r=ones(bombnumber,1); 
for k=l:1:1 
1 = 1 ; 

while i<=bombnumber 

X (i)=randn*sigmax; 
y(1)=randn*sigmay; 
r (i) =sqrt (x(i)''2+y(i)''2); 

1 = 1 + 1 ; 

end 

end 

CEP=median(abs(r)) 

3. Code for Table 4 

Code to calculate the Pk with aiming error (oto), 
ballistic dispersion (rtr), warhead radius 
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clear 

%Ocasion-to-ocasion sigma 
sigmaoto=30; 

%Round-to-round sigma 
sigmartr=10; 

%Warhead Leathal Range 
warheadLR=40; 

%Num of Bombs Per Occasion 
numbpoc=5; 
occ=10; 

kVector=ones(occ,1); 

for i=l:occ %seperate occasions to be simulated 
%xl and yl are the OCO errors 
xl=randn*sigmaoto; 
yl=randn*sigmaoto; 

j=i; 

while j<=numbpoc 

%x2 and y2 are the RTR errors 
x2=randn*sigmartr; 
y2=randn*sigmartr; 

%xt and yt are total error vectors (one col. per bomb) 
xt(j)=xl+x2; 
yt(j)=yl+y2; 

rMiss(j)=(xt(j)"2+yt(j)"2)^.5; 

j=j+i; 

end 


%If rMiss is inside tgt dimensions 
%then a hit has occured 
if min(rMiss)<=warheadLR 
kVector(i)=1; 

else 

kVector(i)=0; 

end 

end 

%Calc Pk from kVector 

Pk=sum(kVector)/length(kVector) 

B. CHAPTER III CODE 


1. Code for Table 5-Table 9 

Investigation into similarity of OTO and RTR as split sigmas vs. one 
total sigma 
clear 

%Ocasion-to-ocasion sigma 
sigmaoto=210 
%Round-to-round sigma 
sigmartr=100 
numoc=1000000; 

%Number of bombs per occasions 
numbpoc=7; 
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for q=l:1:10 

%Split sigmas calculation 

clear xl x2 x3 yl y2 y3 xt yt rl r2 

xt=ones(numoc*numbpoc,1); 

yt=ones(numoc*numbpoc,1); 

r2=ones(numoc*numbpoc,1); 

x3=ones(numoc*numbpoc, 1) ; 

y3=ones(numoc*numbpoc, 1) ; 

rl=ones(numoc*numbpoc, 1) ; 

bombnumber=0; 

1 = 1 ; 

while i<=numoc 

xl=randn*sigmaoto; 

yl=randn*sigmaoto; 

j=i; 

while j<=numbpoc 

bombnumber=bombnumber+l; 
x2=randn*sigmartr; 
y2=randn*sigmartr; 
xt(bombnumber)=xl+x2 ; 
yt(bombnumber)=yl+y2; 

r2 (bombnumber) = (xt (bombnumber) ''2+yt (bombnumber) ^2) 

j=j+i; 

end 

1 = 1 + 1 ; 

end 

cepTwoSigmas(q)=median(r2) ; 

DEPTwoSigmas(q)=median(abs(xt)) ; 

REPTwoSigmas(q)=median(abs(yt)) ; 

%Total sigma calculation 

sigmat=(sigmaoto^2+sigmartr^2)^.5; 

k=l; 

while k<=bombnumber 

x3 (k)=randn*sigmat; 

y3 (k)=randn*sigmat; 

rl(k)=(x3(k)^2+y3(k)^2)".5; 

k=k+l; 

end 

cepOneSigma(q)=median(rl) ; 

DEPOneSigma(q)=median(abs(x3)); 

REPOneSigma(q)=median(abs(y3)); 

end 

cep2avg=mean(cepTwoSigmas) 
ceplavg=mean(cepOneSigma) 
deltaCep=ceplavg-cep2avg 
CEPeql4=l.1774*sigmat 

REP2avg=mean(REPTwoSigmas) 

REPlavg=mean(REPOneSigma) 
deltaREP=REPlavg-REP2avg 
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REPeql2=0.6745*sigmat 


DEP2avg=mean(DEPTwoSigmas) 
DEPlavg=mean(DEPOneSigma) 
deltaDEP=DEPlavg-DEP2avg 
DEPeql3=0.6745*sigmat 

percenterrorCEP=deltaCep/cep2avg 
percenterrorREP=deltaREP/REP2avg 
percenterrorDEP=deltaDEP/DEP2avg 


2. Code for Table 10-Table 11 

Code to calculate the Pk with aiming error (oto), 
ballistic dispersion (rtr), warhead radius 

clear 

%Ocasion-to-ocasion sigma 
sigmaoto=30; 

%Round-to-round sigma 
sigmartr=10; 

%Warhead Leathal Range 
warheadLR=40; 

%Num of Bombs Per Occasion 
numbpoc=5; 
occ=10; 

kVector=ones(occ,1); 

for i=l:occ %seperate occasions to be simulated 
%xl and yl are the OCO errors 
xl=randn*sigmaoto; 
yl=randn*sigmaoto; 

j=i; 

while j<=numbpoc 

%x2 and y2 are the RTR errors 
x2=randn*sigmartr; 
y2=randn*sigmartr; 

%xt and yt are total error vectors (one col. per bomb) 
xt ( j)=xl+x2; 
yt(j)=yl+y2; 

rMiss ( j) = (xt ( j) ''2+yt ( j) ''2) . 5; 

j=j+i; 

end 


%If rMiss is inside tgt dimensions 
%then a hit has occured 
if min(rMiss)<=warheadLR 
kVector(i)=1; 

else 

kVector(i)=0; 

end 


end 
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%Calc Pk from kVector 

Pk=sum(kVector)/length(kVector) 


%THIS CODE IS THE INCORRECT METHOD FOR WEAPON EFFECTIVENESS% 
clear 

%Ocasion-to-ocasion sigma 
sigmaoto=50 
%Round-to-round sigma 
sigmartr=5 

%Warhead Leathal Range 

warheadLR=50 

%number of occasions 

numoc=100000 

%RSS of errors 

sigmat=sqrt (sigmaoto''2 + sigmartr''2) 

%desired pk 
pkdesired=0.7 

%number of rounds to power up (this is the incorrect method!!1) 
ndesired=10 
for q=l:10 

rMiss=zeros(numoc,1); 

for i=l: numoc%seperate weapons to be simulated 
%xl and yl are the OCO errors 
xt=randn*sigmat; 
yt=randn*sigmat; 
rMiss(i)=(xt^2+yt^2)^.5; 
i=i+l; 

end 

%Calc Pk from kVector 

Pk(q)=sum(rMiss<=warheadLR)/length(rMiss) ; 

end 

Pkavg=mean(Pk) 

%now calc required number of bombs to achieve pk=0.7 solve power-up 
formula 

n=log(1-pkdesired)/log(1-Pkavg) 

PkpoweredUP=l-(1-Pkavg)^ndesired 

%THIS CODE IS THE INCORRECT METHOD FOR WEAPON EFFECTIVENESS% 


C. CHAPTER IV CODE 


1. Code for Table 12-Table 25 

% Code to calculate the Pk with oto, rtr, tgtDim (square and circle) 
clear 

%Ocasion-to-ocasion sigma 
sigmaoto=164.469 
%Round-to-round sigma 
sigmartr=164.469 
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%Target Square Dimension Conversion from warhead lethal raduis 

warheadLR=150 

tgtX=warheadLR*pi''. 5 

tgtY=warheadLR*pi^.5 

%Calculations for half tgt side length used to determine a hit 
tgthalfX=tgtX/2 
tgthalfY=tgtY/2 

%One Salvo or number of occasions 
numoc=200000 
for numbpoc=7 
for q=l:10 

clear xt yt kVectorsquare kVectorcircle rMiss 

xt=zeros(numbpoc,1); 

yt=zeros(numbpoc,1); 

rMiss=zeros(numbpoc,1); 

kVectorsquare=zeros(numoc,1); 

kVectorcircle=zeros(numoc, 1); 

for i=l:numoc %seperate occasions to be simulated 
%xl and yl are the OCO errors 
xl=randn*sigmaoto; 
yl=randn*sigmaoto; 

j=i; 

while j<=numbpoc 

%x2 and y2 are the RTR errors 
x2=randn*sigmartr; 
y2=randn*sigmartr; 

%xt and yt are total error vectors (one col. per bomb) 
xt ( j)=xl+x2; 
yt(j)=yl+y2; 

rMiss ( j) = (xt ( j) ''2+yt ( j) ''2) . 5; 

j=j+i; 

end 

%If absolute value xt and yt are inside tgt dimensions 
%then a hit has occured 

kVectorsquare(i)=and(min(abs(xt))<=tgthalfX, min(abs(yt))<=tgthalfY) ; 
kVectorcircle(i)=min(rMiss)<=warheadLR; 
end 

%Calc Pk from kVector 

Pksquare=sum(kVectorsquare)/length(kVectorsquare) ; 

Pkcircle=sum(kVectorcircle)/length(kVectorcircle) ; 

%store results from the above 
PkVectorsquare(q)=Pksquare; 

PkVectorcircle(q)=Pkcircle; 

end 

%average a number of runs to calculate final expected Pk 
Pkavgsquare(numbpoc,1)=mean(PkVectorsquare); 

Pkavgcircle(numbpoc,1)=mean(PkVectorcircle) ; 
end 

Pkavgsquare 

Pkavgcircle 


clear 

c=l % correction factor between .9 and 1 ??? 
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Rt=200 

sigmax=50 %rtr dispersion 
sigmau=150 %oco aiming error 
P=zeros(20, 1) ; 

P2=zeros(20,1); 
for n=l:20 
%eqn 20-11 

for i=l:n 

Pn=(-1) ^(i + 1) ^factorial(n)/(factorial(i)*factorial (n- 
i) ) * (1/i) * ( (Rt''2/sigmax''2) / (2*c+Rt''2/sigmax''2) ) (i-1) * (1-exp (- 

1*(Rt^2/sigmax^2)/((2*c+Rt^2/sigmax^2)/(c*i)+2*sigmau^2/sigmax^2))) ; 
P (n) =P (n) +Pn; 

end 

%eqn 20-13 

for i=l:n 

Pn2=(-1) '' (i + 1) *factorial (n) / (factorial (i) *factorial (n- 
i))*(Rt"2/(sigmax"2*(2+Rt"2/sigmax^2)))"(!- 

1) * ( (Rt''2/sigmax''2 ) / (2+Rt''2/sigmax''2 + i*2*sigmau''2/sigmax''2 ) ) ; 

P2(n)=P2(n)+Pn2; 

end 

end 

P 

P2 


2. Code for Table 26 

% Code to calculate the Pk with oto, rtr, tgtDim (square and circle) 
clear 

%Ocasion-to-ocasion sigma 

sigmaoto=210 

%Round-to-round sigma 

sigmartr=100 

%calculate RSS sigma 

sigmarss= (sigmaoto''2 + sigmartr''2 ) ''. 5 

%Target Square Dimension Conversion from warhead lethal raduis 

warheadLR=220 

tgtX=warheadLR*pi'' . 5 

tgtY=warheadLR*pi^.5 

%Calculations for half tgt side length used to determine a hit 
tgthalfX=tgtX/2 
tgthalfY=tgtY/2 

%One Salvo or number of occasions 
for numbpoc=l:10 
for q=l:10 

clear xl x2 yl y2 xt yt kVectorsquare kVectorcircle rMiss 
for 1=1:20000 %seperate occasions to be simulated 

j=i; 

while j<=numbpoc 

%xt and yt are total error vectors (one col. per bomb) 
xt ( j)=randn*sigmarss; 
yt ( j)=randn*sigmarss; 
rMiss ( j) = (xt ( j) ''2+yt ( j) ''2) . 5; 

j=j+i; 
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end 

%If absolute value xt and yt are inside tgt dimensions 
%then a hit has occured 

kVectorsquare(i)=and(min(abs(xt))<=tgthalfX,min(abs(yt))<=tgthalfY) ; 
kVectorcircle(i)=min(rMiss)<=warheadLR; 
end 

%Calc Pk from kVector 

Pksquare=sum(kVectorsquare)/length(kVectorsquare); 
Pkcircle=sum(kVectorcircle)/length(kVectorcircle) ; 

%store results from the above 
PkVectorsquare(q)=Pksquare; 

PkVectorcircle(q)=Pkcircle; 

end 

%average a number of runs to calculate final expected Pk 
Pkavgsquare(numbpoc,1)=mean(PkVectorsquare); 

Pkavgcircle(numbpoc,1)=mean(PkVectorcircle); 
end 

Pkavgsquare 

Pkavgcircle 


D. CHAPTER V CODE 


1. Double Normal Approximation 

% This program generates DN distributed data with noise and then fits a 
double normal 

% distribution in order to recover the original distribution parameters 

%-curve fitting test program- 

clear 

% First create the data. 
t=-100:2:100; 

t=t(:); % To make t a column vector 
wl=. 7 

sl=30;s2=5; 

Data=l/2* (wl* (1+erf (t/(sl*2''.5))) + (l- 

wl)*(1+erf(t/(s2*2^.5))))+0.02*rand(size(t)); 

figure (1) 

plot(t,Data, 'o ') ; %plot data to be fitted 
ylim([0, 1]) 
hold on 

% Now call FMINSEARCH. 

Starting=[0.1 5 40]; 
options=optimset( 'Display' , 'iter' ); 

Estimates=fminsearch(@DN_myfit,Starting,options,t,Data); 

% To check the fit 
wlel=Estimates(1) 
slel=Estimates(2) 
s2el=Estimates(3) 
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Data_predl=l/2*(wlel*(1+erf(t/(slel*2".5)))+(l- 

wlel) * (1+erf (t/(s2el*2''. 5) ) ) ) ; % function fitted to data 

plot(t,Data_predl, 'r' ); 

xlabel('Miss distance (ft) ') 

ylabel( 'CDF' ) 

grid on 

legend (' Dataset CDF', 'Curve Fitting Results ',' Location ',' SouthEast' ) 
hold off 

% Weapon effectiveness calculation using extracted values and weight 

Wet=20 

Let=20 

numbpoc=l 

for q=l:10 
clear xl kVector 

for 1=1:20000 %seperate occasions to be simulated 
%xl and yl are the OCO errors 
xl=randn*slel; 

%If absolute value xl is inside Wet dimension 
%then a hit has occured 
if abs(xl)<=Let/2 
kVector(i)=1; 

else 

kVector(i)=0; 

end 

end 

%Calc Pk from kVector 

Pkl=sum(kVector)/length(kVector) ; 

%store results from the above 
PkVector1(q)=Pkl; 

end 

%average a number of runs to calculate final expected Pk 
SSPDxl=mean(PkVectorl); 

SSPDxlw=SSPDxl*wlel; 

for q=l:10 
clear x2 kVector 

for 1=1:20000 %seperate occasions to be simulated 
%xl and yl are the OCO errors 
x2=randn*s2el; 

%If absolute value xl is inside Wet dimension 
%then a hit has occured 
if abs(x2)<=Let/2 
kVector(i)=1; 

else 

kVector(i)=0; 

end 

end 

%Calc Pk from kVector 

Pk2=sum(kVector)/length(kVector); 

%store results from the above 
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PkVector2(q)=Pk2; 


end 

%average a number of runs to calculate final expected Pk 
SSPDx2=mean(PkVector2); 

SSPDx2w=SSPDx2*(1-wlel) ; 

SSPDxTotal=SSPDxlw+SSPDx2w 

%Setup Parameters 
mu=0; 
s=slel; 
a=-Wet/2; 
b=Wet/2; 

%Calculate the probability a sample lies within the given range [a:b] 
syms x; 

normpdf=l/ (s*sqrt (2*pi) ) *exp (- ( (x-mu) . ''2) / (2*s''2) ) ; 

Prob=int(normpdf,a,b); 

SSPDlanalytical=double(Prob) ; 

SSPDlanalyticalw=SSPDlanalytical*wlel; 

mu=0; 
s=s2el; 
a=-Wet/2; 
b=Wet/2; 

%Calculate the probability a sample lies within the given range [a:b] 
syms x; 

normpdf=l/ (s*sqrt (2*pi) ) *exp (- ( (x-mu) . ''2) / (2*s''2) ) ; 

Prob=int(normpdf, a, b) ; 

SSPD2analytical=double(Prob) ; 

SSPD2analyticalw=SSPD2analytical*(1-wlel) ; 
SSPDxtotalanalytical=SSPDlanalyticalw+SSPD2analyticalw 


function sse=myfit(params,Input,Actual_Output) 
wlf=params(1); 
slf=params(2); 
s2f=params(3); 

Fitted_Curve=l/2* (wlf* (1+erf (Input/ (slf*2'' .5) ) ) + (1- 
wlf) *(1+erf(Input/(s2f*2^.5)))); 
Error_Vector=Fitted_Curve - Actual_Output; 

% When curvefitting, a typical quantity to 
% minimize is the sum of squares error 
sse=sum (Error_Vector .''2 ) ; 

% You could also write sse as 
% sse=Error_Vector(:)'*Error_Vector(:); 
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Double Rayleigh Approximation 


% This program generates DR distributed data with noise and then fits a 
double normal 

% distribution in order to recover the original distribution parameters 

%-curve fitting test program- 

clear 

% First create the data. 
t=0:2:100; 

t=t(:); % To make t a column vector 
wl=. 7 

sl=30;s2=5; 

Data=l-wl*exp(-t.*t/(2*sl*sl))-(1-wl)*exp(- 
t.*t/(2*s2*s2))+0.02*rand(size(t)); 
figure(1) 

plot(t, Data, 'o ') ; %plot data to be fitted 
ylim([0,1]) 
hold on 

% Now call FMINSEARCH. 

Starting=[0.1 5 40]; 
options=optimset( 'Display' , 'iter' ); 

Estimates=fminsearch(@DR_myfit,Starting,options,t,Data); 

% To check the fit 
wlel=Estimates(1) 
slel=Estimates(2) 
s2el=Estimates(3) 

Data_predl=l-wlel*exp(-t.*t/(2*slel*slel))-(1-wlel)*exp(- 

t. *t/(2*s2el*s2el)); % function fitted to data 

plot(t,Data_predl, 'r' ) ; 

xlabel('Miss distance (ft) ') 

ylabel( 'CDF' ) 

grid on 

legend (' Dataset CDF', 'Curve Fitting Results ',' Location ',' SouthEast' ) 
hold off 

% Weapon effectiveness calculation using extracted values and weight 

Wet=20 

Let=20 

LethalRadius= ( (Wet*Let) /pi) ''0.5 

numbpoc=l 

for q=l:10 
clear xl kVector 

for 1=1:20000 %seperate occasions to be simulated 
%xl and yl are the OCO errors 
xl=randn*slel; 
yl=randn*slel; 
rl=sqrt (xl''2+yl''2) ; 

%If absolute value xl is inside Wet dimension 
%then a hit has occured 
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if rl<=LethalRadius 
kVector(i)=1; 

else 

kVector(i)=0; 

end 

end 

%Calc Pk from kVector 

Pkl=sum(kVector)/length(kVector); 

%store results from the above 
PkVectorl(q)=Pkl; 

end 

%average a number of runs to calculate final expected Pk 
SSPDrl=mean(PkVectorl) ; 

SSPDrlw=SSPDrl*wlel; 

for q=l:10 
clear x2 kVector 

for 1=1:20000 %seperate occasions to be simulated 
%xl and yl are the OCO errors 
x2=randn*s2el; 
y2=randn*s2el; 
r2=sqrt (x2''2+y2''2) ; 

%If absolute value xl is inside Wet dimension 
%then a hit has occured 
if r2<=LethalRadius 
kVector(1)=1; 

else 

kVector(1)=0; 

end 

end 

%Calc Pk from kVector 

Pk2=sum(kVector)/length(kVector) ; 

%store results from the above 
PkVector2(q)=Pk2; 

end 

%average a number of runs to calculate final expected Pk 
SSPDr2=mean(PkVector2); 

SSPDr2w=SSPDr2*(1-wlel) ; 

SSPDrTotal=SSPDrlw+SSPDr2w 

%Setup Parameters 
mu=0; 
s=slel; 
a=0; 

b=LethalRadius; 

%Calculate the probability a sample lies within the given range [a:b] 
syms r 

raypdf= (r/s''2) *exp ( (-r''2) / (2*s''2) ) ; 

Prob=int(raypdf,a,b); 

SSPDlanalytical=double(Prob) ; 

SSPDlanalytlealw=SSPDlanalytleal*wlel; 
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mu=0; 
s=s2el; 
a=0; 

b=LethalRadius; 

%Calculate the probability a sample lies within the given range [a:b] 
syms r; 

raypdf=(r/s^2)*exp((-r^2)/(2*s^2)); 

Prob=int(raypdf,a,b); 

SSPD2analytical=double(Prob) ; 

SSPD2analyticalw=SSPD2analytical*(1-wlel) ; 
SSPDrtotalanalytical=SSPDlanalyticalw+SSPD2analyticalw 


function sse=myfit(params,Input,Actual_Output) 
wlf=params(1); 
slf=params(2); 
s2f=params(3); 

Fitted_Curve=l-wlf.*exp(-Input.*Input/(2*slf*slf))-(l-wlf) .*exp (- 
Input.*Input/(2*s2f*s2f)) ; 

Error_Vector=Fitted_Curve - Actual_Output; 

% When curvefitting, a typical quantity to 
% minimize is the sum of squares error 
sse=sum(Error_Vector.^2); 

% You could also write sse as 
% sse=Error_Vector(:)'*Error_Vector(:); 
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